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FOREWORD 


RECRAM  develcjpnent  was  motivated  by  the  B-IB  project  office  at 
Aeronautical  Systems  Division,  Wri^t-Patterson  AFB,  CH.  Ihe  majority  of  code 
develcpnent  work  was  completed  vdiile  the  authors  were  assigned  to  the  System 
Survivability  Branch  (ENSSS)  in  the  Systems  Analysis  Division  (ENSS)  of  the 
Directorate  for  Systems  Engineering  (ASD/ENS) .  We  thank  Mr.  James  Surifces 
(ENSS)  and  Mr.  Hu^  Griffis  (ENSSS)  for  their  valuable  guidance  and  si^port 
during  the  develcpnent  and  application  of  this  model. 
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EXECUTIVE  SUMMARY 


RESmiAL  DUST  AND  RADIATIC*!  MDDEL  (REERAM) 

MDDEL  DESCKCFTrCX^ 

REaZRAM  is  a  ocnpiter  model  that  simulates  air±)ome  nuclear  cloud 
environments.  Develcpient  was  motivated  by  requirements  for  survivability 
assessments  of  advanced  strategic  edrcraft.  Survivability  assessments  need 
estimates  of  airborne  particle  sizes,  mass  and  radioactivity  levels. 

REEEWM  starts  with  a  stabilized  cloud  of  falling  trace  particles  and 
transports  them  with  spatially  varying  winds  in  a  layered  atmosfAiere.  At  a 
specified  time  after  burst,  REIsaM  determines  the  particle  sizes  that  remain 
aloft;  their  positions  define  the  cloud  axis.  The  figure  below  shows  an 
example  of  a  cloud  axis  at  10  hours  after  a  one  megaton  nuclear  burst. 


Cloud  Axis  From  a  1  Megaton  Weapon,  10  Hours  After  Burst 

REERAM  si;5)eriiiposes  three-dimensional  normal  distributicxi  functions  along 
the  cloud  axis,  vising  particle  positions  as  modes.  The  resulting  distribution 
function  assembly  forms  a  three-dimimTsiQnal  cloud  model  that  can  be  sampled 
along  hypothetical  edrcraft  fli^t  paths. 

REEE^  uses  methods  that  are  supported  by  eitpiriccil  data  v^erever 
possible.  The  variable  wind  fields  are  modeled  with  polyncmial  fits  to 
observed  global  winds.  The  particle  size  distribution  is  based  on  data  from 
atmospheric  nuclear  tests,  as  are  the  cloud  heists,  total  dust  mass  and 
radioactivity  levels.  Validation  ’exercises  included  comparisons  with  Mount 
St.  Helens  ash  cloud  movement  and  ashfall  data.  In  addition,  the  three- 
dimensional  cloud  mass  distributions  were  coaipared  to  aircraft  sample  data 
from  a  hi^  esplosive  dust  cloud  and  to  calculations  from  a  hydrocode. 

The  model  has  been  vised  to  estimte  particle  mass  and  radioactivity  levels 
downwind  of  single  and  multiple  nuclear  bursts.  REEE^AM  is  currently  used  to 
assist  with  the  establishment  of  nuclear  criteria  for  advanced  systems. 
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RECRRM  MDDEL 


RESIDUAL  DUST  AND  RADIATiai  MDDEL  (REDRAM) 


I.  INTRODUCnCW 

REDRAM  is  a  FORTRAN  cxaipiter  code  that  simulates  the  transport  and 
dispersion  of  particles  from  an  atmospheric  nuclear  burst.  The  code  was 
developed  to  estimate  the  dust  particle  sizes,  mass  concentrations  and  radio¬ 
activity  levels  that  are  airborne  long  after  the  burst.  Model  develcpnnent  was 
motivated  by  a  requirement  to  estimate  the  dust  and  radiation  burden  to 
advanced  strategic  aircraft  instrumentation  and  crew.  In  the  event  of  a 
nuclear  war,  aircraft  may  have  to  operate  in  environments  containing  the 
radioactive  d^ris  lofted  weapons  from  intercontinental  and  suhmarine- 
laimched  ballistic  missiles.  Escaping  bombers  and  tankers,  national  command 
authority  aircraft,  and  command,  control  and  communication  aircraft  are 
exaitples  of  systems  that  may  potentially  be  affected  the  nuclear  clouds. 

A  variety  of  coitputer  codes  have  been  developed  to  esqjlore  late  time 
airborne  dust  and  radiation  environments  (13;  19;  22).  The  models  have  been 
vised  to  estimate  environments  with  varying  degrees  of  fidelity,  depending  v^ion 
their  intended  use.  Fast  running  codes  sacrifice  resolution  for  speed  of 
operation;  they  are  intended  for  studies  that  cannot  tolerate  the  long  run 
times  and  tedious  interpretation  required  ty  the  larger,  more  detailed 
methods.  However,  all  conputer  models  of  nuclear  clouds  suffer  from  the  fact 
that  there  is  no  atmospheric  nuclear  test  data  that  will  svpport  their 
coDtplete  and  vinambigious  validation.  ^ 

The  R&D  Associates  (RDA)  dust  and  radiation  environment  model  (19)  was 
developed  for  Headquarters,  Defense  Nuclear  Agency  (ENA) .  RDA  based  their 
model  on  as  much  reliable,  enpiriceil  data  as  possible.  RDA  used  the  model  for 
a  general  assessment  of  environment  severity  following  a  hypothetical  large 
scale  nuclear  exchange.  ENA  and  RDA  cocperated  with  REEI^'s  develcpment  by 
providing  valuable  references  and  constructive  comments  regarding  the 
technical  areas  vhere  inprovements  would  be  most  beneficial. 

RECRAM  is  philosophically  similar  to  the  RDA  model  in  that  conputational 
methods  were  based  on  enpirical  correlations  of  observed  test  data  vherever 
reliable  data  was  available.  REERAM  significantly  differs  from  the  RDA  model 
in  at  least  two  respects:  RECRAM  uses  a  continuously  varying  qpectrally 
deccmposed  wind  field  for  inproved  lateral  transport  calculations,  and  REERAM 
development  has  stressed  validation  of  tran^iort  methods  with  as  much  relevant 
test  data  as  available  in  order  to  develop  confidence  in  simulation  fidelity. 

This  report  describes  the  REERAM  code,  and  its  submodels  of  nuclear  envi¬ 
ronments.  Efforts  to  validate  the  wind-driven  transport  simulation  are  summa¬ 
rized.  Results  of  a  one  megaton  cloud  transport  calculation  illustrate 
typical  cloud  translation,  shear  and  airborne  mass/radiation  gradients. 
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n.  MDDEL  DESC3?IPnON 


General 

REODRAM  starts  with  a  stabilized  cloud  of  trace  particles,  tran^rts  them 
downwind  until  they  fall  to  the  ground,  then,  backing  vp  in  t±ne,  reconstructs 
a  three  dimensional  model  of  dust  and  radiation  environments  at  a  xoser- 
specif ied  time.  Of  course,  the  specified  time  must  precede  the  time  at  vhich 
the  smallest  trace  particle  lands. 

The  stabilized  doud  is  a  gravity-sorted  array  of  user-defined  trace 
particle  sizes.  Different  particle  sizes  start  falling  frm  different 
altitudes,  d^)ending  on  weapon  yield.  Starting  heists  are  based  on  ancilyses 
of-  stabilized  wafer  heists  in  the  D^iartment-  of  Defense  Land  Fallout 
Interpretive  Code  (DELFIC)  (10;  16) .  One  particle  per  size  is  used  in  REDRAM 
to  represent  the  mode  of  a  distribution  of  same-size  particles  in  the 
vertical.  Each  falling  trace  particle  is  tracked  from  its  starting  hei^t, 
throu^  a  variable  wind  field,  to  the  ground.  Intermediate  particle  positions 
are  recorded,  and  used  to  interpolate  back  to  find  the  position  of  each 
particle  at  a  user-defined  freeze  frame  time.  The  connection  of  all  airborne 
particles  at  freeze  frame  time i; is  the  cloud's  nearly  vertical  freeze  frame 
axis,  skewed  by  the  spatially  varying  wind  field  throu^  vMch  the  trace 
particles  fell. 

The  particle  positions  are  used  as  modal  locations  for  ej^anding  three- 
dimensional  normal  distributions.  The  distributions  from  all  particles  are 
superinposed  at  atmospheric  layer  heists  vp  and  down  the  freeze  frame  axis. 
The  result  is  an  axially  skewed  cloud  of  overlapping  density  functions  in 
three  dimensions.  The  normal  density  functions  ^proximate  the  way  cloud  mass 
density  and  radiation  vary  with  location  about  the  cloud's  axis.  The  space 
around  the  cloud  is  divided  into  horizontal  medies,  vertically  stacked  with 
user-defined  resolution. 

ComDutational  Me^ 

RECRAM  divides  the  atmo^here  into  a  user-defined  number  of  vertical 
layers.  The  ipper  and  lower  boundaries  of  layers  are  called  levels,  so  a  ten 
layer  model  has  eleven  levels.  Atmospheric  prcperties  (density,  kinematic 
viscosity)  are  ocnputed  at  each  level.  Winds  are  defined  on  a  different  set 
of  vertical  levels;  they  are  specified  by  the  National  Meteorological  Center 
at  twelve  different  ischaric  heists  above  sea  level.  Winds  are  ccstputed  at 
the  iscbaric  heists  that  are  above  and  belcw  particle  starting  heists,  then 
the  winds  are  interpolated  to  the  particle  heists.  Particles  start  from 
size-dependent  heists  in  the  initial  cloud,  and  fall  throui^  atmospheric 
layers  until  they  reach  ground.  Thus,  initially,  there  are  three  simultaneous 
schemes  for  representing  altitude-varying  parameters:  1)  the  coiputational 
mesh  of  evenly  divided  layers,  2)  the  iscbaric  wind  levels  and  3)  initial 
heists  of  particles  in  the  stabilized  cloud.  Figure  1  illustrates.  REDRAM 
interpolates  winds  and  particles  onto  the  conputational  mesh  in  the  process  of 
conputing  the  pxositipns  of  falling  trace  particles. 
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ATM)SPHER£  CLOUD  SPECTRAL 

LEVELS  AND  PARTICLE  HEIGHTS 

LAYERS  DIA.  (wa)  (MILLIBARS) 


PARnCIE  HEIGHTS 
ARE  FOR  A  ONE 
MEGAim  dDOD 


Figure  1.  Ccmpitatlonal  Me^,  Wind  Levels  and  Initial  Particle  Heists 


it  Data 


The  viser  sv:5plies  information  alx3ut  the  burst,  the  particle  size 
frequency  functions  and  model  paramet^  settings.  Specifically,  Table  1 
identifies  the  data  necessary  for  a  REES^AM  ccmputer  run.  Ir^wts  are  contained 
in  two  separate  units,  one  for  the  spectral  coefficients  and  one  for  all  other 
information. 


TABLE  1 


REDRAM  INTOT  DATA  REQUIREMENTS 

Weapon  Yield 
Burst  Hei^t 
Trace  Particle  Sizes 
Activily  Fraction  in  each  Particle  Size 
Mass  Fraction  in  each  Particle  Size 
Fission  Fraction 
Total  Mass  Lofted 
Particle  Mass  Density 
Latitude  and  Longitude  of  Burst  Point 
Numbers  of  Standard  Deviations  that  Define 
Limits  of  Three  Dimensional  Me^ 

Lateral  Spatial  Resolution 
Wind  Spectral  Coefficients 
Freeze  Frame  Time 
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stabilized  Clcjud  Heicflit 


Starting  heists  of  particles  in  the  initial,  stabilized  cloud  are  set  by 
fits  to  the  average  heists  of  wafers  in  DEUTC  code  oirtput.  Reference  10 
explains  how  the  following  yield-d^)endent  fits  were  cbtained.  Ctonstants  are 
in  Table  2. 

Zp  =  “S  Bd  +  I 

In  S  =  SO  +  SI  Iny  +  S2  (Iny)^  +  S3  (lny)3  +  S4  (Iny)^ 

In  I  =  10  +  II  Iry  +  12  (lny)2  +  I3  (Iny)^  +  I4  (Iny)^ 

vAiere:  Zp  =  particle  hei^t  (meters) 

Bd  =  particle  diameter  (micrcroeters) 
y  =  weapon  yield  (kilotons) 

Sn  and  In  =  n-th  constant 

lny.=  natural  logarithm  of  y 


TABLE  2 


BOLYNCMIAL  CONSTANTS  FOR  EARTICIE  STARTING  HEIGHTS 
n  Sn  In 


0 

1 

2 

3 

4 


1.574 

-0.1197 

0.03636 

0.00410 

0.0001965 


7.889 

0.3400 

0.001226 

0.005236 

0.000417 


The  correlations  give  particle  heists  that  are  gravity  sorted;  that  is, 
larger  particles  start  ftcm  Icwer  altitudes  than  smaller  particles.  The 
particles  are  assumed  to  be  spheres  falling  thzouc^  a  U.S.  Standard 
Atmosphere. 

Atmosphere 

Atmospheric  properties  are  needed  to  determine  each  particle's  hei^t- 
dependent  fall  speeds  and  residence  tiroes  in  the  model's  vertical  layers. 
REI3RAM  uses  the  U.S.  Standard  Atmosphere  generating  equations  (14)  to  ooatrpute 
atmospheric  state  variables  needed  for  terminal  velocity  ccarputation.  The 
code  determines  U.S.  standard  air  mass  density  and  kinematic  viscosity  at  the 
top  and  bottom  of  each  layer.  The  U.S.  Standard  Atmosphere  is  used  to 
approximate  conditions  at  mid-latitudes  in  the  northern  hemisphere.  Sensi¬ 
tivity  studies  with  atmospheric  variables  from  trcpical  latitudes  show  that 
particle  fall  speeds  are  not  very  sensitive  to  atinospheric  prop)erties  (10) . 
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Particle  Fall  Speeds 


The  t-P>rminaT  velcxiity  of  each  trace  particle  size  is  ccaipited  at  the 
and  bottcm  of  each  atinospheric  layer,  and  the  average  is  vised  to  determine 
that  particle's  residence  time  in  the  layer. 

tr  =  dz/Vz 

vhere;  tr  =  residence  time 

Vz  =  terminal  velocity 
dz  =  layer  thickness 


REERAM  uses  the  methods  of  Davies  and  McDonald  (6;  12)  to  ccropute  the 
terminal  velocities  of  Uniform  density  spheres.  Ihe  Davies-IfecDonald  methods 
are  caramon  to  several  nuclear  fallout  codes  (e.g.,  references  3  and  16) . 
calculations  are  valid  for  spherical  particles.  If  the  itajority  of  particles 
are  not  ^eres,  then  results  must  be  interpreted  as  equivalent  sphere 
transport  and  transfer  functions  would  be  required  to  derive  particle  size- 
spiecific  masses  and  radiation  doses. 


Wind  Vector  Components 

Ihe  lateral  translation  of  falling  particles  in  each  layer  is  performed 
by  multiplying  each  particle's  residence  time  by  the  wind  speed  at  the  center 
of  the  layer  at  the  particle's  latitude  and  longitude.  It  is  assumed  that 
wind  speed  and  direction  are  constant  in  a  model  layer  only  for  the  span  of 
time  that  a  particle  resides  there.  ^ 

Wind  vector  ocnponaits  are  ocmputed  with  a  set  of  polynanials  for  each 
conponent  at  each  of  twelve  ischaric  heists.  Ihe  pxjlyncanials  ar&  fits  to 
glc±al  gridded  wind  data.  Ihe  fits  are  spherical  harmonics  that  describe  the 
variation  of  each  wind  ccaotponent  as  a  function  of  latitude  and  longitude  (9) . 
Each  polynomial  is  a  finite  sum  of  spectral  terms  that  describe  the  vector 
conponent.  Spectral  deccarposition  accurately  r^licates  the  wave-lite 
variations  of  wind  speed  around  the  earth.  A  general  form  of  the  spectral 
wind  expansion  follows: 


vhere: 


U,V 


A 

J 


n 


P 

1 


wind  vector  ccnponents 

ocarplex  spectral  coefficients  for  U,V 

1  vhen  SI  =  0,  otherwise  A  =  2 

truncation  limit 

latitudinal  index 

longitudinal  index 

associated  Legendre  polynomial 

latitude 

longitude 

square  root  of  -1 
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The  spectral  coefficients  are  ccnputed  external  to  REEEffiM  with  software 
provided  by  the  Air  Force  Geophysics  laboratory  (APGL) .  Monthly  most  probable 
wind  data  is  gridded  (2  1/2 ‘  x  2  1/2' )  and  fit  with  ^erical  harmonics.  The 
fits  yield  the  spectral  coefficients.  The  above  equations,  with  spectral 
coefficients,  permit  the  calculation  of  wind  vector  ccnaponent  magnitude  at  any 
latitude  and  longitude. 

This  method  is  based  on  the  National  Meteorological  Center's  twelve  layer 
30  wave  spectral  model  of  the  atmosphere.  There  are  961  (31  x  31)  spectral 
coefficients  required  for  each  of  the  two  wind  vector  coirponents  at  each  of 
the  12  isobaric  levels. 

Monthly  most  prciable  winds  were  provided  ky  the  USAF  Environmental 
Technical  i^lications  Center  (USAFEIAC) .  Winds  were  chosen  by  analyzing 
seven  years  of  global  wind  data,  and  selecting  specific  days  with  the  surface 
features  and  vpper  air  flow  patterns  most  typical  of  each  month.  USAFEIAC 
provided  gridded  data  for  monthly  most  typical  days  and  for  one  day  before  and 
after  each  typical  day.  USAFEIAC 's  data  was  further  resolved  by  interpolation 
to  obtain  a  2  1/2'  x  2  1/2'  wind  field.  Those  gridded  winds  were  then  used  as 
iiput  to  the  APGL  software  to  obtain  spectral  coefficients.  Winds  were  then 
conputed  from  the  coefficients  at  grid  points  and  comtpared  to  USAFEIAC  data  to 
veri^  the  coefficient  set. 

Wind  Shear 

The  lateral  growth  of  a  nuclear  cloud  is  estimated  with  equations  that 
require  values  of  wind  shear.  Shear  is  defined  in  REDEEM  as  follows: 


sx,  ^  = 


vhere:  sx,sy 
Au,  Av 


Az 

tr 

tfall 

N 


^ear  components  (x  or  y) 

change  in  wind  speed  (u  or  v  component) 
over  a  vertical  distance,  Az 

layer  thioikness 

particle  residence  time  in  layer  j 
total  fall  time  of  particle 

total  number  of  layers  throun^  vhich  a  particle  falls 


Shear  is  a  residence  time  wei^ted  average  of  shears  that  the  fcilling  cloud 
experiences  in  each  of  the  model  layers.  Shear  will  be  positive,  expanding 
the  cloud  as  time  progresses. 
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Cloud  Lateral  Growth 


As  the  cloud  falls,  different  particle  sizes  fall  with  different  speeds, 
so  they  are  translated  by  winds  for  different  lengths  of  time.  It  is 
possible,  and  typical,  for  the  tcp  and  bottcm  of  a  clcxid  to  be  in  different 
locations  at  a  given  time.  Ihe  cloud  breads  laterally,  first  with  initial 
toroidal  motion,  then,  with  winds  that  ^ear  tiie  cloud.  Smeill  scale  atmo¬ 
spheric  motions  dilute  the  cloud  cind  diffuse  ttie  particles  away  from  the 
cloud's  axis.  Ihe  cloud's  lateral  growth  is  a  very  coaiplicated  process  to 
simulate  precisely,  and  attempts  to  do  so  have  met  with  varied  degrees  of 
success  (1) .  EEDRAM  uses  the  lateral  growth  model  that  was  developed  for  the 
WSBG-10  fallout  code  (20)  and  corrected  by  Bridgnan  in  the  AFTT  fallout  code 
(2;  3) .  Ihe  lateral  growth  model  includes  ej^licit  d^>endenoe  on  weapon  yield 
and  wind  shear. 

Ihe  AFTT  code  is  a  constant  wind  fallout  model  that  employs  a  semi- 
empirical  ej^ression  for  crosswind  cloud  dispersion.  Ihe  method  was  adapted 
for  PEDES'S  variable  wind  treatment  because  the  method  includes  the  most 
influential  causes  of  cloud  growth,  and  because  the  method  was  developed  to 
fit  atmospheric  test  cloud  data.  Ihe  AFTT  equations  estimate  the  value  of  a 
standard  deviation  of  an  assumed  normal  probability  density  function  oriented 
in  the  crosswind  direction.  In  RECRAM,  the  same  equations  are  xased  to  give 
directional  standard  deviations  with  directional  shear  values  computed  from 
^jectral  winds. 

Cloud  Vertical  Growth 

In  addition  to  the  cloud  ^reading  caused  by  variable  fall  rates  of 
different  sized  particles,  cloud  thioikness  is  affected  by  toroidal  flew  and 
^trainment/dif fusion  processes.  However,  because  vertical  velocil^  wind 
^ears  are  generally  so  small  compared  to  horizontal  velocity  shears,  toroidal 
growth  and  particle  fall  dominate  the  process  of  vertical  growth.  Particle 
fall  is  ej^licitly  treated  with  trace  particle  tracking.  Trace  particle 
positions  are  assumed  to  be  the  modes  of  normal  distributions  aligned  with  the 
cloud's  freeze  frame  axis.  Standard  deviations  of  those  modes  are  measures  of 
cloud  vertical  growth,  again  based  on  fits  to  DELFIC  cloud  calculations  (5) . 
RECRAM  models  a  cloud's  vertical  standard  deviation  as  25%  of  the  visible 
cloud  thickness.  Vertical  standard  deviations  vary  with  weapon  yield  and 
particle  size  in  the  initial  cloud.  Ihe  vertical  Gaussian  distribution 
associated  with  each  particle  size  falls  at  a  rate  that  is  determined  by  the 
trace  particle  size,  overlapping  in  time  and  space  with  distributions  from 
other  particle  sizes. 

Clrmd  Mac;s  loading 

An  atmospheric  nuclear  burst  can  raise  a  large  mass  of  ddDris  if  the 
burst  occurs  close  to  the  ground.  Reference  15  contains  a  d^iction  of  anhy¬ 
drous  cloud  mass  loading  vs  burst  hei^t  normalized  by  weapon  yield.  Ihe  plot 
is  reproduced  in  Figure  2. 
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Figure  2.  Cloud  Mass  Loading 

The  cloud  mass  loading  curve  was  digitized  and  fit  with  five  quadratics. 
In  (M/Y)  =  ao  +  ai(x)  +  33  (x)  2 

v^iere:  IJ/Y  -  negatons  (MT)  of  d^ris/megaton  of  weapon  yield 

X  =  In  (SHOE) 

SFTOB  =  scaled  hei^t  of  burst  (feet/MrV^) 

In  =  natural  logarithm 

Table  3  shows  the  values  of  quadratic  constants. 


TABLE  3 

Quadratic  Constants  for  Mass  Loading  Equations 


X 

ao 

ai 

^2 

X  <  2.2215 

-  1.31 

-.048 

0 

2.2215  OC  <  7.9537 

-  2.984 

1.042 

-0.1707 

7.9537  <X  <  8.3337 

-82.637 

22.258 

-1.580 

8.3337  <X.  <  8.8289 

-325.684 

84.066 

-5.498 

8.8289  ^ 

-18141.701 

4148.995 

-237.351 
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Mass  loading  is  an  iirprecisely  known  quantity.  Data  is  sparse,  contro¬ 
versial  and  of  questionable  relevance  to  bursts  over  soil  types  that  are  dif¬ 
ferent  frcm  desert  alluvium.  However,  several  of  the  mass  loading  estimtes 
were  derived  by  independent  techniques,  agreeing  on  the  approxhnate  1/3 
megaton  of  dust  per  megaton  of  yield  for  surface  bursts.  Uhoertainty 
increases  with  burst  hei^t. 

Particle  size  Distributions 

A  recent  survey  of  particle  size  distributicai  (PSD)  functions  illiKtrated 
the  wide  variation  in  PSDs  that  can  be  produced  by  nuclear  bursts  (22) .  ESDs 
vary  with  the  same  features  that  mass  loadings  do:  soil  type,  water  coitent, 
rock  strength  and  burst  hei^t  to  name  a  few.  The  REESWM  user  can  specify  ary 
PSD  by  supplying  trace  particle  sizes  and  the  mass  and  activity  density 
function  values  associated  with  them.  However,  to  be  consistent  with  RECRAM 's 
development  philosophy  of  using  enpirically  based  information,  the  DELFIC- 
default  number  size  spectrum  is  currently  used  for  code  development  and 
verification. 

DELFIC-default  is  a  peuticle  distribution  derived  frcm  measurements  of 
particle  sizes  that  fell  frcm  the  cloud  produced  by  a  lew  yield  atmos^eric 
nuclear  test  in  Nevada.  The  DELFIC-default  activity-size  distribution  is 
accurately  replicated  as  the  wei^ted  sum  of  two  log-normal  density  functions. 


vhere:  a(r)=  fraction  of  particles  with  radius  r 
r  =  particle  radius 

fv  =  fraction  of  the  population  with  radioactivity 
distributed  within  the  particle's  volume 

1-fv  =  fraction  of  the  population  witii  radioactivity 
distributed  on  the  particle's  surface 

^  =  log  slope  of  distribution 

<^3/  ^^2  ~  medians  of  the  volume  and  surface  distributions 

Values  of  the  activity  fraction  and  mass  feaction  are  inputs  to  REERAM,  not 
hard-wired,  so  that  other  PSDs  can  be  studied. 
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Output  Retxjrts 


EECRAM  output  echoes  the  irput  data  except  the  spectral  coefficients, 
vdiich  are  too  numerous  to  list  each  time  that  REI»AM  cperat^.  Particle 
trajectories  are  summarized  by  printing  particle  position  coordinates  at  each 
level  in  the  atmosphere.  locations  of  points  on  the  freeze  fr^  cloud  suds 
are  printed  along  with  shears  and  conversions  to  earth  radians  in  a  spherical 
coordinate  system.  The  output  listing  includes  a  level  by  level  expansion  of 
the  three  dimensional  normal  distribution  function  values  for  each  particle 
size.  The  expansions  are  r^rted  after  PECRAM  r^jositions  the  freeze  frame 
cloud  axis  onto  user-defined  grid  rode  locations.  The  e3?)ansions  are  printed 
first  for  the  northeast  quadrant  of  a  cloud,  then  values  are  r^rted  again 
for  the  fully  ejqanded  cloud.  Full  eiqpansion  uses  symmetry  properties  of  the 
distribution  functions  to  reflect  values  into  the  other  three  cloud  quadrants. 
Figure  3  illustrates  the  process  of  translating  the  particles,  freezing  their 
positions  at  a  user-defined  time,  oonputing  density  function  values  in  the 
northeast  quadrant,  then  reflecting  the  northeast  quadrant's  nuDt±iers  into  the 
other  quadrants. 


Figure  3.  REEE?AM  Qaraputation  Sequence 
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III.  ninr.T.  ENVIPONMEWr  CAIOJIAnONS 

Using  the  previously  described  methods  of  ocnputing  relevant  cloud 
bdiavior,  REESIAM  determines  tmit  time  reference  dose  rate  and  particulate  mass 
concentrations  in  each  cell.  Following  is  a  general  description  of  the 
model's  approach  to  calculating  cell  environments. 

Dose  Rate 

REERAM  reports  values  of  unit  time  reference  dose  rate  and  particle  mass 
concentrations  at  discrete  points  in  the  atmosphere.  Ihose  points  represmt 
the  centers  of  three  dimensional  cells  that  are  one  layer  thickness  hi^ 
(meters)  and  one  increment  of  lateral  mesh  (radians)  in  each  horizontal 
direction.  Dose  rate  is  skyshine  from  particles  within  the  cell  volume,  so 
pall  dimensions  should  be  larger  than  the  mean  free  paths  of  gamma  photons 
emitted  by  radioactive  elements  in  the  particles.  Figure  4  illustrates  the 
gecmetry. 


Figure  4.  CJell  -  Receiver  Gecmetry 
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The  dose  rate  to  a  receiver  located  at  cell  center  is  the  sum  of  each 
particle  size's  activity,  converted  to  dose  rate,  attenuated  by  intervening 
air  and  multiplied  by  REI»AM's  determination  of  how  much  of  the  three  dimen¬ 
sional  density  function  is  in  the  cell. 


VolUMC. 


e _ 

4-trs^  di'^olwMe. 


vhere:  = 

NER  = 
Ar  = 

Ykt  = 
f  f  = 
F  = 
SNC  = 
= 

<h'^>  = 

M.I.  = 


s 


unit  time  reference  dose  rate 

number  of  different  particle  sizes  in  cell 

finite  particle  radius  interval 

weapon  yield 

fission  fraction 

REI»fiM  calculation  of  three  dimensional  density  function 

source  normalization  constant 

mass  attenuation  coefficient  for  receiver 

average  energy  of  photons  emitted  by  radioactive  particles 

total  mass  attenuation  coefficient  for  air 

mass  integral  of  air  between  source  of  radiation  and  receiver  ‘ 

path  length  traveled  by  photon 


REDRAM  assumes  530  roegacuries  of  gamma  radioactivity  per  kiloton  of  fission 
(8) .  Ihe  volxmne  integral  is  solved  analytically  by  assuming  that  the  over- 
vhelming  majority  of  radiation  is  coning  fron  within  the  cell,  so  the  integra¬ 
tion  can  be  performed  within  an  infinite  sphere.  Typical  values  of  attenua¬ 
tion  coefficients  were  used  for  air  and  receiver  materials  (7) . 
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since  the  dose  rate  values  are  for  one  hour  after  burst  (unit  time  refer¬ 
ence) ,  that  Vcdue  most  be  integrated  from  receiver  arrival  time  to  receiver 
departure  time  to  obtain  dose  to  receiver. 

-i.ii 

ji 

vhere:  =  dose  (rads)  to  receiver 

ta  =  arrival  time 

td  =  d^arture  time 

Assuming  that  the  enseaiible  of  fission  products  decay  with  a  net 
dependence,  the  above  integral  can  be  solved  analytically  for  total  dose  to 

the  receiver.  This  st^  is  not  performed  in  REEEffiM,  since  arrival  and 

departure  times  are  scenario  d^iendent. 

Particulate  Mass 

Ihe  total  particulate  mass  in  a  cell  is  the  sum  of  the  masses  contributed 
by  each  particle  size.  REES^AM's  three  dimensional  spatial  distribution 
function  gives  the  fraction  of  each  particle  size  that  is  contained  in  each 
cell.  Multiplying  that  fracticai  by  the  fraction  of  mass  contained  in  each 
particle  size,  then  by  the  total  cloud  mass  yields  total  mass  contributed  to  a 
cell  ty  each  particle  size. 


M(r)  =  (TM)  (F)  m(r) 


Vihere:  M(r)  =  total  mass  ccartributed  to  a  cell  hy  particle  size  r 

IM  =  total  cloud  mass  lofted 

F  =  RECRAM  ccilculation  of  three  dimensioned,  density  function 

m(r)  =  fraction  of  IM  contained  in  trace  particle  size  r 

The  values  of  m(r)  for  each  trace  particle  size  are  irpjts  to  REDRftM.  The 
total  cloud  mass,  TM,  comes  from  the  quadratic  fits  to  the  mass  loading  curve. 
Figure  2. 
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IV.  SaMPIE  CALjCUIATICNS 


REDRAM  was  vised  to  cxjipute  the  cloud  particle  mass  and  radioactivity 
levels  associated  with  a  one  megaton  nuclear  burst.  Following  is  a  list  of 
initial  conditions  for  these  sairple  calculations. 


Weapon  Yield:  1  megaton 
Burst  Hei^t:  0  (surface  burst) 

Mass  Fraction  Per  Particle  Size  =0.03 

Fission  Fraction:  0.5 

Total  Mass  Lofted:  0.33  megatons 

Particle  Mass  Density  :  2.6  grams  per  cubic  centimeter 
Latitude  and  Longitude  of  Burst  Point:  40N,  lOOW 
Lateral  Spatial  Resolution  =  0.05  degrees 
Wind  Spectral  Coefficients:  16  July  1981,  1200Z 

Figures  5  and  6  shew  cloud  mass  density  contours  at  8  and  10  hours  after 
burst.  These  sequential  slices  of  cloud  clearly  illustrate  the  effects  of 
variable  winds  translating  and  shearing  the  cloud.  Figures  7  and  8  illustrate 
the  cloud's  axis  at  the  same  two  times.  Variable  wind  effects  are  again 
evident,  as  is  the  cloud's  general  motion  downwind.  Trace  particles  have  also 
fallen  to  lower  altitudes. 

REDRAM  was  also  used  to  ooitipute  the  airborne  dust  mass  densities  produced 
ty  50  simultaneous  one  megaton  bursts.  The  calculations  simulated  a  hypothe¬ 
tical  attaede  on  a  missile  field  in  the  north  central  Ifriited  States.  Single 
burst  calculations  were  linearly  superimposed  to  estimate  the  multiple  burst 
environments.  Reference  4  describes  this  REDRAM  ^plication  for  reentry 
vehicle  fratricide  studies. 
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LRTITUOE  (DEGREES)  LATITUDE  (DEGREES) 

38.73S  39.382  «.029  10.676  11.323  11.970  12.617  13.261  38.73S  39.382  10.029  10.676  11.323  11.970  12.617  13.261 
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Figure  5.  Cloud  Density  Contours,  8  hours 


4000  METERS  (  10  HRS  )  —  MASS  DENSltY  (MILLIGRAMS/Mxm3) 
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Figure  6.  Cloud  Density  Contours,  10  Hours 
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Figure  7.  Cloud  Axis  at  8  Hours 


Figure  8.  Cloud  Axis  at  10  Hours 
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V. 


VAT.imTICN  SUMMRRY 


Since  nuclear  cloud  siniulation  requires  so  mich  estimation  and  modeling 
with  so  little  unambiguous  data,  a  significant  portion  of  development 

work  was  spent  on  validation  exercises.  Frcm  the  beginning,  code  development 
philosophy  was  to  vise  the  best  available  empirical  models.  This  section 
briefly  summarizes  the  results  of  three  different  validation  zeroises. 
REDEEM  was  intended  for  late  time  environment  calculations  to  svpport 
survivability  assessments  of  advanced  strategic  eiircraft.  Accordingly,  much 
emphasis  was  placed  on  development  of  the  c^ability  to  transport  the  clouds 
realistically,  with  variable  winds,  for  long  times.  Validation  studies  of  the 
new  spectral  wind  transport  methocis  were  performed  with  two  recent,  hi^ily 
reliable  data  sources  and  a  hydrocode  calculation. 

First,  the  spectral  wind  transport  method  was  used  to  ccsipute  the 
time-varying  position  of  the  adi  cloud  lofted  by  the  18  Ifey  1980  eruption  of 
the  Mount  St.  Helens  volcano.  That  ccaiparison  confirmed  that  the  basic 
approach  of  moving  a  gravity-sorted  particle  stack  in  a  Lagrangian  reference 
frame  can  replicate  the  motion  of  a  real  cloud  in  the  atmosphere  with 
surprisingly  good  fidelity  (9;  11) . 

Second,  REDKAM  was  adapted  to  cxmpute  the  time-varying  position  and  con¬ 
tent  of  the  dust  cloud  created  by  the  raRECT  COURSE  hi^  ejplosive  (HE)  test. 
DIRECT  0CXIE?SE  was  a  600  ton  elevated  spherical  charge  that  was  detonated  at 
White  Sands  Missile  Range,  New  Mexico  on  26  October  1983.  The  cloud  was 
tracked  and  dust  measurements  were  taken  by  an  instrumented  aircraft  at  a 
variety  of  altitudes.  The  REI®fiM  code,  with  the  HE  cloud  hei^t  and  lateral 
cloud  spread  substituted  for  the  nuclear  cloud  correlation  was  exercised  to 
compute  cloud  position  and  airborne  mass  densities  (9;  17)  at  times  that 
corresponded  to  sampling  eiircraft  passes.  A  particle  size  distribution  was 
developed  frcm  ground  core  sample  data  and  used  in  REDE^  to  compute  specific 
mass  density  values.  Trace  particle  sizes  were  selected  to  correspond  to  size 
bins  in  the  aircraft  sampling  instrumentation,  and  results  were  path 
integrated  for  a  direct  ccnparison  with  aircraft  instrumentation  output.  The 
DIRECr  axJRSE  cloud  position  was  accurately  ccnputed  by  REIS^  with  spectral 
coefficients  for  26  October  1983  svpplied  by  the  National  Meteorological 
Center.  The  cloud  mass  densities  also  ccnpared  favorably  with  the  aircraft 
sample  data.  REE^AM  histograms  of  mass  fraction  per  size  range  as  a  function 
of  particle  size  agreed  with  aircraft  data  best  \^en  samples  were  taken  at 
hi^  enou^  altitudes  to  preclude  orographic  effects  on  the  wind  field. 
Reference  17  contains  a  complete  sumimciry  of  the  REI®AM  veilidation  exercise 
with  DIRECT  COURSE  data.  Generally,  that  effort  confirmed  the  quality  of 
REDRAM 's  spectral  wind  transport  method  and  produced  quantitative  evidence 
that  airborne  dust  mass  levels  can  be  estimated  accurately. 
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Third,  RECRAM  results  were  cxatpared  to  calculations  made  with  the 
Defense  Nuclear  Agency's  DICE  (Dirt  Inplicit  Ctode  Eulerian)  hydrocode  (21). 
DICE  is  the  agency's  beiseline  code  for  simulating  details  of  cloud  formation 
and  growth.  RECRAM  results  were  in  good  agreement  with  the  more  scphisticated 
and  detailed  DICE  data.  For  example.  Figure  9  illustrates  the  similar 
contours  of  lOE-8  grams  per  cubic  centimeter  particle  mass  density  computed  by 
REERAM  and  DICE.  A  more  detailed  description  of  the  REERAM-DICE  comparison  is 
contained  in  Reference  18. 


1  MT  Surface  Burst  at  10  Minutes 


Ronge,  kft 


Figure  9.  REIRAM  -  DICE  Cloud  Ccaiparison 
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